#ifndef JACOBIN_HH
#define JACOBIN_HH

extern void Jacobin(const double* q, const double* a, const double* d, double* J_p, double* J_r);
extern void Mmt(const double* A, double* B);
extern void Mmm(const double* A, double* B, double* C);
extern void Mcond(const double* eigval, double* cond);
extern double Dvec2(const double* a, const double* b);
extern void QRdecop(const double* A, double* Q, double* R);
extern void sorteig(double* A);

#endif
